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Abstract: Leguminosae are an important part of terrestrial ecosystems and play a key role in promoting 
soil nutrient cycling and improving soil properties. However, plant composition and species diversity 
change rapidly during the process of succession, the effect of leguminosae on soil physical-chemical and 
biological properties is still unclear. This study investigated the changes in the composition of plant 
community, vegetation characteristics, soil physical-chemical properties, and soil biological properties on 
five former farmlands in China, which had been abandoned for 0, 5, 10, 18, and 30 a. Results showed that, 
with successional time, plant community developed from annual plants to perennial plants, the importance 
of Leguminosae and Asteraceae significantly increased and decreased, respectively, and the importance of 
grass increased and then decreased, having a maximum value after 5 a of abandonment. Plant diversity 
indices increased with successional time, and vegetation coverage and above- and below-ground biomass 
increased significantly with successional time after 5 a of abandonment. Compared with farmland, 30 a of 
abandonment significantly increased soil nutrient content, but total and available phosphorus decreased 
with successional time. Changes in plant community composition and vegetation characteristics not only 
change soil properties and improve soil physical-chemical properties, but also regulate soil biological 
activity, thus affecting soil nutrient cycling. Among these, Leguminosae have the greatest influence on soil 
properties, and their importance values and community composition are significantly correlated with soil 
properties. Therefore, this research provides more scientific guidance for selecting plant species to 
stabilize soil ecosystem of farmland to grassland in the Loess Plateau, China. 
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1 Introduction 


Secondary succession is a key measure to restore degraded ecosystems and plays an important 
role in increasing soil nutrient content and plant diversity (Zeng et al., 2017). During secondary 
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succession, plant productivity and community composition are affected directly (Li et al., 2019; 
Liu et al., 2019). Some studies have shown that plant characteristics (such as plant biomass, 
coverage, evenness, and diversity) increased after secondary succession (Fan et al, 2015; 
Tessema and Belay, 2017), and soil properties may increase with succession due to the increase in 
plant biomass (Wang et al., 2009; Zhao et al., 2019). Some researchers have shown that natural 
regeneration through secondary succession can restore degraded soil properties and maintain soil 
fertility (Li et al., 2019, 2022; Zou et al., 2022; Li et al., 2023). Soil physical-chemical properties 
are affected by decomposition of plant root exudates and litter decomposition matter due to 
vegetation growth (Zhao et al., 2015). In addition to soil microbial biomass, ecological enzyme 
activity is sensitive to changes in soil and plants. Therefore, secondary succession affects soil 
physical-chemical properties, enzyme activity, and microbial biomass (Zhu et al., 2012). 

It has been shown that long-term vegetation succession is a key driver of succession in restored 
areas (van der Heijden et al., 2008; Lucas-Borja et al., 2019; Hao and Chu, 2021), and that 
vegetation changes inevitably affect soil nutrient (Reich, 2005; Peichl et al., 2012). Moreover, 
fast-growing vegetation can affect soil nutrient storage through litter decomposition and root 
exudates (Berg, 2014; Lozano et al., 2014). van der Putten et al. (2013) found that plants could 
influence soil microbial communities and drive changes in soil physical-chemical properties. In 
addition, plants can improve microclimate conditions and increase soil microbial activity, thus 
promoting soil nutrient cycling and accumulation (Wickander et al., 2021; Bao et al., 2022). In 
contrast, soil extracellular enzymes are rate-limiting steps of microbial metabolism, affecting 
microbial decomposition of carbon and nutrients, and their activity is also positively correlated 
with organic carbon content, affecting soil chemical properties (Chabrerie et al., 2003), but the 
pattern of change in enzyme activity during long-term natural processes of secondary succession 
is unclear. And different plant species have different effects on soil because each plant has its own 
unique biological characteristics that may affect nutrient return and decomposition of organic 
carbon differently (Fanin and Bertrand, 2016). For example, Leguminosae typically contain high 
levels of nitrogen, while C4 plants contain high levels of carbon and can significantly increase 
nitrogen and carbon contents of soil (Sivaram et al., 2018; Yao et al., 2021). Some studies have 
shown that rapid decomposition of litter of nitrogen-rich plants can increase soil nutrient content 
(Horodecki and Jagodzinski, 2017). In addition to interspecific differences, it has also been found 
that effects of the same plant functional group on soil properties vary depending on successional 
time (Hu et al., 2016). Thus, the effects of changes in vegetation on soil properties are variable 
and uncertain during secondary succession (Zhang et al., 2016; Zhang et al., 2022). 

There are fewer studies on the effects of Leguminosae on soil properties (Fterich et al., 2014), 
and the potential mechanisms of their plant communities on soil properties need to be determined 
to help evaluate the soil nutrient balance and predict the sustainability of restoration (Song et al., 
2022; Liao et al., 2023). The Loess Plateau of China, which covers an area of about 64x10^ kny?, 
has characterized by severe soil erosion and excessive desertification (Zheng, 2006; Zhang et al., 
20182). The Chinese government has undertaken a series of ecological restoration projects such 
as returning farmland to forest and grassland (Fu et al., 2000), which is an important measure to 
increase vegetation coverage and prevent fragmentation and degradation of ecosystem (Lozano et 
al, 2014; Zhang et al, 2016). After the converting of farmland into forest and grassland, 
secondary succession without human disturbance, biodiversity, and soil ecosystem functions have 
been significantly improved. Meanwhile, soil physical-chemical properties (Qiu et al., 2022), 
microbial community dynamics (Liu et al., 2022), soil enzyme activities (Chen et al., 2022), plant 
community composition (Zhu et al, 2021a), and plant diversity (Zhang, 2005) have been 
significantly changed. However, there are few researches on the variability of Leguminosae on 
soil physical-chemical properties and biological properties. In this study, we investigated the 
dynamic plant community composition, plant diversity, soil physical-chemical properties, and soil 
biological properties using different years of secondary successional abandoned cropland and 
farmland as the subjects. We hypothesized that plant community composition and plant diversity 
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would affect soil physical-chemical properties and soil biological properties to different degrees 
with succession time. In addition, we further hypothesized that Leguminosae are the dominant 
species in the plant community and have a greater influence on soil properties during secondary 
succession. Therefore, we focused on: (1) effects of succession time on plant community 
composition, species diversity, soil physical-chemical properties, and soil biological properties; (2) 
effects of plant community characteristics and community composition on soil properties; and (3) 
effects of Leguminosae on soil properties. 


2 Materials and methods 


2.1 Study area 


A field experiment was conducted in the Zhifanggou watershed of Ansai District (36°42'-36°47'N, 
109?13'—109?16'E), Shaanxi Province, northwestern China, which is in the center of the Loess 
Plateau. The soil in this area is classified by the Food and Agriculture Organization (FAO) as 
extremely erodible Calcaric Cambisol that developed from wind-blown loess deposits. Due to 
long-term wind-blown deposits and soil erosion, the catchment (8.73 km?) has the typical 
landforms of the loess hilly and gully landscape with slope gradients varying from 0? to 65?, and 
altitudes ranging from 1000 to 1300 m a.s.l. (Wang et al., 2009). Climate in this region is 
classified as temperate with an average annual precipitation of 510 mm, of which approximately 
70% falls between June and September (Zhang et al., 2016). The average frost-free period in this 
catchment is 7150 d and the annual average temperature is 8.8?C, with the maximum temperature 
of 36.8°C in August and the minimum temperature of —23.6°C in January. The Zhifanggou 
watershed serves as a template for remedying soil erosion and ecological environment. In 1973, 
the Chinese government began to carry out the integrated management of soil and water 
conservation and vegetation restoration. One of the most widely used methods to prevent soil 
nutrient loss is to convert former farmlands into abandoned grassland for natural recovery without 
anthropogenic interference (Zhang et al, 2016). After decades of this artificial restoration, 
vegetation coverage in the area increased significantly, and large-scale herbaceous communities 
formed, which included the species Artemisia sacrorum Ledeb., Lespedeza dahurica Schindler, 
Artemisia scoparia Waldst. et Kit, and Heteropappus altaicus (Willd.) Novopokr. In addition, 
there are still a handful of farmlands in the gentle slope zone where maize, potatoes, and foxtail 
millet have been planted. 


2.2 Experimental design 


This study used the space-for-time substitution method, which has been applied to numerous field 
studies (Walker et al., 2010; Williams et al., 2013; Zhang et al., 2016; Zhang et al., 2018b) and a 
meta-analysis (Zhou et al., 2017) to investigate the effects of natural succession on the above- and 
below-ground ecosystem. Through a comprehensive field survey of the Zhifanggou watershed in 
Ansai District, four plots were selected. Farmlands that had been abandoned 5, 10, 18, and 30 a 
were randomly chosen from areas with similar gradients, slopes, and altitudes. Active farmland 
was used as the control (CK). Detailed information about each succession sequence of abandoned 
farmland and referential farmland is presented in Table S1. All selected sample plots had been 
plowed every spring before abandonment, and were planted mainly with maize and cereals, with 
minimal fertilizer application and irrigation relying mainly on rainfall. There was no human 
disturbance during the abandonment period and vegetation was naturally colonized and growing. 
Each succession sequence of the abandoned farmland was represented by three independent 
replicates and space between any two sites was large enough to exclude spatial dependence («14 
m) for the vast majority of soil and plant variables (Marriott et al., 1997). Three replicated plots 
(30 mx30 m; Fig. 1) were randomly set up in each plot for subsequent investigation and sampling. 
For each variable, average value of the three replicates in the same plot was taken to constitute 
observations. Thus, in total, 15 observations were established (5 succession sequencesx3 
replicates) for each variable. 
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Sampling plot 


Fig. 1 Sampling process of the Zhifanggou watershed in the Loess Plateau, China 


2.3 Vegetation investigation and soil sampling 


Vegetation investigation and soil sampling were conducted in August 2022. Nine replicated 
sample points were selected along an "S" shape in each plot for soil sampling (Fig. 1). After 
removing the litter horizons and biological crusts, nine soil samples of 0-10 cm depth were 
collected from each point using a 5-cm diameter stainless steel auger, then fully homogenized to 
provide one composite sample per plot. These soil samples were immediately sieved through «2 
mm mesh to remove visible litter, animal residue, roots, stones, and debris. Part of each soil 
sample was immediately taken back to the laboratory and stored at 4?C for biological 
characteristics analysis, and other subsamples were air-dried and stored at room temperature for 
physical-chemical analysis. Soil water content (SWC) and bulk density (BD) were obtained 
randomly from six points per plot using a 5-cm diameter steel core sampler (Fig. 1), then dried at 
105?C to constant weight. 

Vegetation characteristics were investigated in situ, and 10 quadrats with an area of 1 m? were 
randomly arranged in each plot to record the name, height, coverage, frequency, and density of 
each species. Plant coverage was estimated at each plot. Both above- and below-ground parts of 
all plants in each 1 mx1 m quadrat were collected by clipping and digging, respectively. After 
washing the roots with tap water, both above- and below-ground parts of all plants were dried at 
75°C to constant weight to calculate above-ground biomass (AB) and below-ground biomass 
(BB). The species importance value (IV), Margalef richness index (M), Shannon-Wiener diversity 
index (H), and Pielou evenness index (E) were calculated (Zhang, 2005; Zhang et al., 2016): 

IV=(relative density-relative frequency+relative coveragetrelative height)/4x 100%, (1) 
where relative density (%), relative frequency (%), relative coverage (%), and relative height (%) 
are defined as the percentages of density, frequency, coverage, and height of a single species to 
the total density, frequency, coverage, and height of all species per plot, respectively. 


M -(S-D/lnN, (2) 
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Ss 


H - onn). ) 


i=l 


S 
E= -5 (piln p;)/nS, (4) 
i=l 
where S is the total number of species in each plot; N is the sum density of all species in each 
plot; and P; is the density proportion of species i. 


2.4 Soil physical-chemical properties analysis 


Soil organic carbon (SOC) concentration was measured using K2Cr2O7 oxidation method (Bao, 
2000). Soil total nitrogen (TN) concentration was determined using a continuous flow analyzer 
(AA3, Bran+Luebbe, Norderstedt, Germany) after micro-Kjeldahl digestion (Zhang et al., 2016). 
Concentrations of ammonium nitrogen (NH4*-N) and nitrate nitrogen (NO3-N) of the soil were 
determined using a continuous flow analyzer (AA3, Brant+Luebbe, Norderstedt, Germany) after 
extractions of fresh soil with 2 M KCI for 18 h (Zhang et al., 2016). Soil total phosphorus (TP) 
concentration was analyzed using Mo-Sb anti-Spectro photograohy method (Bao, 2000) and soil 
available phosphorus (AP) was determined by Olsen method (Olsen et al., 1982). Soil alkaline 
phosphatase (ALP), urease (URE), catalase (CAT), and saccharase (SAC) activities were 
determined using the method described in previous study (Zhang et al., 2018b), which is shown in 
Table S2. Soil microbial biomass carbon (MBC), nitrogen (MBN), and phosphorus (MBP) were 
estimated from fresh soil samples using a chloroform fumigation-extraction method (Bao, 2000; 
Ren et al., 2016a). SWC was determined by oven drying samples at 105°C to constant weight 
(Bao, 2000). Soil pH was measured using a pH meter (Sartorius PB-10, Goettingen, Germany) 
after vibrating soil-water (1:5) suspension for 30 min (Ren et al., 2017; Zhang et al., 2018b). Soil 
BD was determined using soil core method, and obtained by calculating the ratio of soil mass to 
total volume after oven-drying at 105°C to constant weight (Ren et al., 2016b). Soil temperature 
was measured using a right-angle geothermometer. Soil clay was determined by a laser particle 
size analyzer. 


2.5 Soil microbial respiration and metabolic quotient measurement 


Soil microbial respiration was routinely determined using alkali (NaOH) absorption method as 
described by both Hu et al. (1997) and Liu et al. (2010). Briefly, 25 g of fresh soil was weighed 
from each plot and placed evenly in a 500-mL glass flask. The blank treatments (without soil) were 
used as controls. A 25-mL flask with 10 mL 0.5 mol/L NaOH solution was placed in the 500 mL 
glass flask to capture CO» evolved from the soil. There apparatuses were incubated at 25?C in the 
darkness for two weeks, and cover of each flask was opened for half an hour once a day in previous 
week and once every two days in the next week. Soil microbial metabolic quotient (qCO2) was 
calculated as MR/MBC (where MR is microbial respiration) (Wardle and Ghani, 1995). 


2.6 Statistical analyses 


Distribution of all data was checked using Shapiro-Wilk test, and all variables were found to 
follow a normal distribution. One-way analysis of variance (ANOVA) and Duncan's multiple 
range test (P«0.05) were used to assess the effects of successional stage on plant community 
composition, plant species diversity, soil physical-chemical properties, and soil biological activity 
using R v.3.1.3 software (R Core Team, 2015). Correlations between plant characteristics, soil 
physical-chemical properties, and soil biological activity were identified using redundancy 
analysis (RDA) and relative interpretation rate of vegetation characteristics on soil properties. 
Venn diagram analysis of relative interpretation rates of Leguminosae dominant species on soil 
physical-chemical and biological properties was undertaken. The clustering of different samples 
along successional stage by variables of soil and plant characteristics were revealed through 
nonmetric multidimensional scaling (NMDS) method. Both RDA and NMDS were conducted by 
CANOCO v.5.0 software package (Ter Braak and Smilauer, 2002). Analysis of interpretation rate 
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of Leguminosae dominant species with soil physical-chemical and biological properties used 
Venn diagram package at R v.4.2.1. 


3 Results 


3.1 Changes in plant community characteristics 


Vegetation characteristics changed significantly with successional time (Table 1). Plant 
community coverage, AB, and BB increased rapidly at 5 a plot, decreased at 10 a plot, and 
significantly increased again from 10 to 30 a plot (Table 1). Plant community coverage, AB, and 
BB varied from 40.44% to 58.94%, 116.12 to 275.06 g/m?, and 44.58 to 117.54 g/m”, respectively, 
and the minimum and maximum values occurred at 10 and 30 a plots, respectively (Table 1). The 
M, H, and E significantly increased with succession time, ranging from 1.62 to 3.51, 1.13 to 3.65, 
and 0.47 to 0.95, respectively (Table 1). 


Table 1 Plant community characteristics over different successional time 


Years of farmland abandonment 


Index F P 
5a 10a 18a 30a 
Coverage (%) 58.8743.82° 40.44+2.47° 51.59+2.61° 58.94+1.70 30.04 <0.001 
Above-ground biomass (g/m?) 145.9945.97° 116.12+8.04? 228.27416.27°  275.06:8.53* 147.06 «0.001 
Below-ground biomass (g/m?) 52.46+4.53° 44.58+4.08° 91.87£7.15° 117.54+8.21° 90.34 <0.001 
Margalef richness index 1.62+0.08° 2.67+0.27° 2.92+40.27° 3.51+0.09* 48.17 <0.001 
Shannon- Wiener diversity index 1.13+0.074 2.50+0.13° 2.89+0.29° 3.65+0.12* 114.38 <0.001 
Pielou evenness index 0.47+0.034 0.65+40.02° 0.80+0.03° 0.95+0.02" 227.99 <0.001 


Note: Different lowercase letters within the same row indicate significant difference among different successional time at P<0.05 level. 
Mean+SD; n=3. 


Plant community composition was considerably affected with succession time (Table 2), and 
NMDS suggested that sampling plots were clearly gathered into four well-differentiated groups at 
plant community level (Fig. S1). IV of Artemisia capillaris Thunb. significantly decreased with 
succession time, and plant species dominated the community at 5 and 10 a plots and turned into a 
companion species for those plots abandoned for more than 18 a. Although H. altaicus was also 
the dominant species at 5 a plot, its dominance was replaced by Stipa bungeana Trin. that appears 
at 10 a plot and was not present at 30 a plot. IV of L. dahurica was significantly increased during 
the first 18 a of abandonment, was decreased at 30-a plot, and dominated the plant community at 
18-30 a plot. Two perennial species, A. sacrorum and Artemisia giraldii Pamp., appeared over the 
course of succession at 10 a plot, and were dominant species at 18 and 30 a plots, respectively. In 
general, dominant species of plant community appeared to gradually transition from A. capillaris 
and H. altaicus to L. dahurica, A. sacrorum, and A. giraldii during 30 a of abandonment. For plant 
species composition, IV of Compositae plants decreased rapidly during the first 18 a of 
abandonment and then increased at 30 a plot, with the lowest point at 18 a plot (31.10%) and the 
peak at 5 a plot (54.09%). IV of Gramineae plants significantly increased and then decreased with 
succession time, and the peak occurred at 18 a plot (29.43%). IV of Leguminosae plants 
significantly increased with successional time, ranging from 9.86% to 36.49%. 


3.2 Changes in soil physical-chemical properties 


NMDS revealed that soil physical-chemical properties significantly changed with successional 
time (Fig. S2). SOC and soil clay was significantly increased, ranging from 3.53 to 5.70 g/kg and 
18.86% to 20.81% for farmland and 30 a plot, respectively (Fig. 2; Table 3). Compared with 
farmland, SWC, TN, NH4-N, and NO3-N contents were significantly lower at 5 a plot but 
significantly higher with successional time and rising to farmland levels at 10 a plot (Table 3). 
SOC, TN, NH4-N, NO;-N, SWC, and clay contents reached their highest values at 30 a plot, and 
were 61.47%, 45.51%, 13.11%, 65.32%, 36.23%, and 10.34% higher than those at farmland, 
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respectively (Table 3). BD, soil temperature, TP, and AP values were higher at farmland than at 
those plots, and decreased significantly with increasing successional time, ranging from 1.15 to 


Table2 Importance value of plant species over different successional time 


Years of farmland abandonment 


Plant species 5a 10a 18a 30a 
(%) 

Patrinia heterophylla Bunge - 4.56+0.49* 1.89:0.26^ 0.58+0.09° 
Dracocephalum moldavica L. 4.03+0.48* - - - 
Lespedeza dahurica Schindler 9.86+0.114 16.36+0.66° 21.20+0.53° 17.79+0.40° 
Medicago sativa L. - - 3.56+0.44° 6.53+0.57" 
Vicia sepium L. - - 2.3340.51* - 
Thermopsis lanceolata R. Br. - - - 3.07+0.64" 
Glycyrrhiza uralensis Fisch. - - 1.22+0.39° - 
Astragalus melilotoides Pall. : 3.8940.28° 2.98+0.62° 6.11+0.46° 
Gueldenstaedtia verna Boriss. - - - 2.99+0.34° 
Setaria viridis (L.) Beauv. 10.29+0.36* 5.36+0.46° 0.68+0.06° - 

Poa pratensis L. 6.41+0.27° 2.09+0.15° 3.3040.33° - 

Stipa bungeana Trin. - 19.41+0.77° 13.45+0.99° 8.60+0.79° 
Cleistogenes chinensis (Maxim.) Keng - - 2.50+0.58° 0.55+0.36° 
Roegneria kamoji Ohwi - - 4.54+0.50* 0.30+0.11° 
Phragmites australis (Cav.) Trin.ex Steud - - - 0.75+0.18" 
Leymus secalinus (Georgi) Tzvel. - - - 1.41+0.13? 
Bothriochloa ischaemum (Linnaeus) Keng - - 4.96+0.65° 12.74+0.64* 
Viola philippica Cav. - 1.2140.27° - - 
Artemisia capillaris Thunb. 35.02+40.44* 18.31+0.78° 3.68+0.18° 2.07+0.76¢ 
Heteropappus altaicus (Willd.) Novopokr 17.37+0.80" 8.29+40.74° 2.17+0.44° - 
Cirsium setosum (Willd.) MB. 1.70+0.14? - - - 
Artemisia sacrorum Ledeb. - 7.8340.72° 16.49+0.90* 10.11+0.62° 
Ixeris polycephala Cass. - 2.30+0.38° - 1.840.32° 
Saussurea japonica (Thunb.) DC. - - 2.59+0.25* 1.2340.25° 
Artemisia giraldii Pamp. - 2.24+0.21° 5.12+0.55° 18.88+0.93* 
Ixeridium sonchifolium (Maxim.) Shih - 1.59+0.13? - 0.50+0.12° 
Bidens parviflora Willd. - - 1.052:0.09* - 
Salsola collina Pall. 4.87+0.06* 1.24+0.26° - 0.8540.11° 
Chenopodium glaucum L. - 2.79+0.78* - - 
Potentilla discolor Bge. 6.99+0.61° 2.53+0.47° - - 
Duchesnea indica (Andr.) Focke - - 3.37+0.43" - 
Polygala tenuifolia Willd. 1.17+£0.33° - 0.57+0.06° 1.86+0.10° 
Incarvillea sinensis Lam. 2.29+0.20° - 2.34+0.38° 1.26::0.06^ 
Compositae 54.09+40.66" 40.57+£0.21° 31.10+0.384 34.6340.22° 
Gramineae 16.70+0.44¢ 26.86+0.17° 29.43x1.22* 24.35x0.26* 
Leguminosae 9.86+0.114 20.25+0.93° 31.30+1.62° 36.49:0.08* 
Sum of Compositae, Gramineae, and Leguminosae 80.65+0.708 87.67+1.27° 91.82+0.63° 95.48+0.59° 


Note: Different lowercase letters within the same row indicate significant difference among different successional time at P<0.05 level. 


Mean+SD; n=3. - means no value. 
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1.29 g/cm?, 16.44?C to 18.71°C, 0.53 to 0.60 g/kg, and 0.61 to 1.95 mg/kg, respectively (Table 3). 
In addition, soil pH was slightly higher at all plots than at farmland, although there was no 
statistical difference between plots (Table 3). 

RDA showed that soil physical-chemical properties were significantly correlated with 
vegetation characteristics and plant community composition (Fig. 3). For vegetation 
characteristics, the first two axes significantly explained 96.3796 of all considered variables, 
indicating that M, H, E, AB, BB, Leguminosae, Gramineae, Compositae and sum of Leguminosae, 
Gramineae, and Compositae were positively correlated with SOC, TN, NH4'-N, NO3-N, SWC, 
and clay, and negatively correlated with TP, AP, BD, and soil temperature (Fig. 3a). Conversely, 
Compositae plants were negatively correlated with SOC, TN, NH4*-N, NO; -N, SWC, and clay, 
and positively with TP, AP, BD, and soil temperature. For plant community composition, the first 
two axes explained 75.68% of all considered variables, and SWC, AP, NH4*-N, and NO; -N were 
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Fig. 2 Changes of soil nutrient contents over different successional time. (a), soil organic carbon (SOC); (b), 
total nitrogen (TN); (c), total phosphorus (TP); (d), available phosphorus (AP); (e), nitrate nitrogen (NO3 -N); (f), 
ammonium nitrogen (NH4'-N). Different lowercase letters indicate significant difference among different 
successional time at P«0.05 level. Bars are standard errors. 


Table3 Soil physical-chemical properties over different successional time 


Years of farmland abandonment 


Parameter Farmland F P 
5a 10a 18a 30a 
Soil water content (76) 7.2440.18°  6.7140.164 7.82::0.13* 9.28::0.26^ 9.862:0.10* 177.30 | «0.001 
pH 8.49:0.02* — 8.5440.03** 8.62::0.02* 8.49+0.03° 8.56+0.03° 10.72 0.001 
Soil bulk density (g/cm?) 1.29+0.02° 1.26+0.02° 1.21+0.03° 1.18+0.02"° 1.15+0.02° 23.30 <0.001 
Soil temperature (°C) 18.7120.24*^ 18.14+0.13° 17.52+0.30° —16.44+0.264 17.56+0.32° 32.16 — «0.001 
Clay (96) 18.86£0.81° 19.2240.28° 19.52+0.26%  20.19+0.33% — 20.8120.29* 9.21 0.002 
SOC (g/kg) 3.5340.319 3.9440.22°¢ 4.3140.23° 4.98+0.28° 5.70:0.20* 35.28 «0.001 
TN (g/kg) 0.4140.02% — 0.3640.05* 0.45+0.02° 0.51+0.02° 0.59+0.03? 32.47 «0.001 
TP (g/kg) 0.6020.01* — 0.5850.01* . 0.57+0.01% 0.53+0.014 0.55::0.01* 14.43 «0.001 
NH4-N (mg/kg) 9.45+0.44> — 7.564:0.38? 8.66::0.42* 9.82+0.30° 10.68+0.56° 23.13 «0.001 
NO; -N (mg/kg) 3.1240.23° 2.4740.124 3.5940.29° 4.47+0.21° 5.15+0.26° 63.54 «0.001 
AP (mg/kg) 1.95x0.17^ — 1.64+0.07° 1.20+40.16° 0.78+0.074 0.610.074 69.98 — «0.001 


Note: SOC, soil organic carbon; TN, total nitrogen; TP, total phosphorus; AP, available phosphorus. Different lowercase letters within 
the same row indicate significant difference among different successional time at P<0.05 level. MeanzSD; n-3. 
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Fig. 3 Redundancy analysis (RDA) for the relationships among soil physical-chemical properties, vegetation 
characteristics, and plant community composition. (a), relationships between soil physical-chemical properties 
(red arrows) and vegetation characteristics (grey arrows); (b), relationships between soil physical-chemical 
properties (red arrows) and plant community composition (grey arrows). SOC, soil organic carbon; TN, total 
nitrogen; TP, total phosphorus; AP, available phosphorus; SWC, soil water content; BD, bulk density; T, soil 
temperature; M, Margalef richness index; H, Shannon-Wiener diversity index; E, Pielou evenness index; AB, 
above-ground biomass; BB, below-ground biomass; sum of C, G, and L, sum of Compositae, Gramineae, and 
Leguminosae. 


the most influential factors causing the changes in plant composition, contributing to 48.40%, 
47.80%, 46.70%, and 44.20% of total variation (Fig. 3b). 


3.3 Changes in soil biological properties and responses to plant community characteristics 


Four soil enzymatic activities (SAC, URE, ALP, and CAT) responded consistently with 
successional time (Figs. 4a-d and S3). Notably, the activities of four soil enzymes were 
significantly higher at abandoned lands than at farmland, with the increased rates ranging from 
25.46% to 292.19% (SAC), 38.08% to 112.34% (URE), 39.02% to 211.1296 (ALP), and 4.23% to 
21.37% (CAT), respectively. All four soil enzymatic activities increased significantly with 
successional time and reached their maximum values at 30 a plot. Similarly, MBC and MBN 
significantly increased with successional time, with growth rates ranging from 26.33% to 
156.55% and from 17.58% to 116.79% for farmlands and 30 a plot, respectively (Fig. 4e and f). 
Soil microbial respiration increased and peaked at 10 a plot and then sharply decreased with 
successional time (Fig. 4g). Soil respiration at abandoned lands was higher than that at farmland 
(Fig. 4g). Value of qCO» increased significantly at 5 a plot compared with farmland, then 
dramatically decreased with successional time (Fig. 4h). In addition, values of qCO» at 18 and 30 
a plots were significantly lower than those at farmland, 5, and 30 a plots. 

RDA showed a significant correlation between soil biological properties and plant community 
characteristics (Fig. 5). For vegetation characteristics, the first two axes significantly explained 
93.37% of all considered variables (Fig. 5a). In particular, Leguminosae plants were the most 
important factor, contributing 81.60% of total variation. RDA indicated that M, H, E, AB, BB, 
Leguminosae, Gramineae, and total of Leguminosae, Gramineae, and Compositae were positively 
correlated with MBC, MBN, SAC, URE, ALP, and CAT and negatively correlated with qCO» (Fig. 
5a). Conversely, Compositae plants were negatively correlated with MBC, MBN, SAC, URE, ALP, 
and CAT, and positively correlated with qCO» (Fig. 5a) In addition, all plant community 
composition variables significantly explained 78.83% of total soil biological variation (Fig. 5b). 


3.4 Effect of Leguminosae on soil physical-chemical properties and biological characteristics 


RDA showed that Leguminosae had the greatest influence on soil physical-chemical properties, 
with a relative explanation percentage of 75.9096, which was significantly higher than other 
vegetation characteristics (Table 4). Moreover, RDA showed that Leguminosae had the greatest 
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Fig. 4 Changes of soil biological properties over different successional time. Different lowercase letters indicate 
significant difference among different successional time at P«0.05 level. Bars are standard errors. C, carbon; N, 
nitrogen. (a), saccharase; (b), urease; (c), alkiline phosphatase; (d), catalase; (e), microbial biomass C (MBC); (f), 
microbial biomass N; (g), soil microbial respiration; (h), metabolic quotient. 


influence on soil biological properties, with a relative explanation percentage of 81.6096, which 
was significantly higher than other vegetation characteristics (Table 5). Among plant community 
compositions, Leguminosae were the most dominant community, with the most dominant species 
Medicago sativa L. making the greatest relative contribution to soil physical-chemical and 
biological properties with 10.60% and 36.90%, respectively (Fig. 6). In summary, Leguminosae 
plants can be regarded as the most important plant factors associated with soil physical-chemical 
and biological properties, but not dominant species and species diversity. 
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Fig. 5 Redundancy analysis (RDA) for relationship among soil biological properties, vegetation characteristics, 
and plant community composition. (a), relationships between soil biological properties (red arrows) and 
vegetation characteristics (grey arrows); (b), relationships between soil biological properties (red arrows) and 
plant community composition (grey arrows); MBC, soil microbial biomass carbon; MBN, soil microbial biomass 
nitrogen; ALP, alkaline phosphatase; URE, urease; CAT, catalase; SAC, saccharase; MR, microbial respiration; 
qCO2, metabolic quotient; M, Margalef richness index; H, Shannon-Wiener diversity index; E, Pielou evenness 
index; AB, above-ground biomass; BB, below-ground biomass; sum of C, G, and L, sum of importance values of 
Compositae, Gramineae, and Leguminosae plants. 


Table4 Effects of vegetation characteristics on soil physical-chemical properties over successional time 


Parameter Contribution (96) F P 
Leguminosae 75.90 27.4 0.002 
Compositae 66.10 17.6 0.002 
Gramineae 37.90 5.8 0.016 
Above-ground biomass 59.40 13.4 0.004 
Below-ground biomass 59.10 13.3 0.002 
Margalef richness index 53.30 10.6 0.002 
Shannon-Wiener diversity index 58.30 12.9 0.002 
Pielou evenness index 56.20 11.8 0.002 


Table 5 Effects of vegetation characteristics on soil biological properties over successional time 


Parameter Contribution (96) F P 
Leguminosae 81.60 31.6 0.002 
Compositae 59.00 13.6 0.006 
Gramineae 27.30 3.6 0.052 
Above-ground biomass 73.70 25.8 0.004 
Below-ground biomass 77.70 31.8 0.002 
Margalef richness index 45.10 7.9 0.020 
Shannon-Wiener diversity index 47.60 8.4 0.004 
Pielou evenness index 46.60 8.4 0.010 


4 Discussion 


4.1 Plant characteristics over successional time 
A. capillaris, an annual species belonging to Compositae, was a pioneer species that first colonized 
the bare farmland and quickly became the dominant species at 5 a plot (Table 2). The 
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Fig. 6 Variance partitioning analysis of dominant species of Leguminosae on soil physical-chemical properties 
(a) and soil biological properties (b) 


strong resource competitiveness and rapid environmental adaptability led to rapid growth of A. 
capillaris and further inhibited the growth of other species (Zhang et al., 2016), resulting in low 
H, M, and E at 5 a plot (Table 1). This study showed that H, M, and E significantly increased with 
successional time until 30 a plot (Table 1), and more species had appeared in 18 and 30 a plots, 
such as S. bungeana, A. sacrorum, and Astragalus melilotoides Pall. (Table 2). Similar results 
were found in previous studies, with high species diversity and richness occurring during middle 
successional stage of the Chinese Loess Plateau (25 a of abandonment) (Wang et al., 2009; Zhang 
et al., 2016; Sun et al., 2017). The colonization by plants improved soil conditions and provided 
more living space for more species during secondary succession (Zhu et al., 2021b; Wang et al., 
2022), and the improvement in soil condition could be demonstrated by increased SWC and clay 
content and decreased BD and pH (Table 3). At the same time, many species of different ecotypes 
interacted with each other to reach a stable state in these improved community environments (Liu 
et al., 2023; Yang et al., 2023), thereby causing an increase in species diversity and richness (Sun 
et al., 2022; Gu et al., 2023; Kong et al., 2023). In addition, the results showed that IV of 
Leguminosae and Gramineae plants increased with successional time (Table 2), indicating that 
these families are more adaptable to the arid environment of the Loess Plateau than others (Zhang 
et al., 2021; Liu et al., 2022). For example, Leguminosae have deeper roots to obtain water and 
nutrients (Hu et al., 2016). Morever, M, H, and E were found to be positively correlated with 
Leguminosae and Gramineae plants (Figs. 3, 5, and S2). It should be noted that strong C and N 
fixation capacity of Gramineae and Leguminosae plants can increase soil C and N input to 
improve availability (Hu et al., 2016; Tian et al., 2016), thereby creating suitable conditions for 
plant colonization to increase species diversity and richness. In contrast, Gramineae and 
Leguminosae plants can coexist with other species due to complementarity, facilitation, and niche 
differentiation (HilleRisLambers et al., 2004; Wu et al., 2017), resulting in increased species 
diversity and richness. However, Compositae plants can inhibit the growth of other species due to 
allelopathic effects (Balicevié et al., 2016), leading to them being negatively correlated with M, H, 
and E (Figs. 3 and 5). Above all, these results validate the hypothesis that plant community 
composition and species diversity significantly changed with successional time. 


4.2 Responses of soil physical-chemical and biological properties to plant community 
characteristics 


Soil physical properties are sensitive to alterations in plant characteristics (e.g., population 
composition, species diversity, and biomass) (Ile et al., 2021; Cotrufo et al., 2022). The results 
indicated that dynamics of plant diversity (M, H, and E) and coverage corresponded to the 
accumulation of above- and below-ground biomass (Table 1), with both being related to SWC, 
BD, soil temperature, and clay (Fig. 3a). This result showed that secondary succession of plant 
communities can alter soil physical properties (Jiao et al, 2011). First, the increase in 
above-ground biomass implies an improvement in vegetation and litter coverage (O'Halloran et 
al., 2013; DeBerry and Atkinson, 2014), thereby reducing solar radiation on the soil surface (Zhao 
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et al., 2015) and ultimately leading to an increase in SWC and a decrease in soil temperature with 
successional time (Table 3). Secondly, community composition changed from annual herbs to 
perennial Compositae, Leguminosae, and Gramineae plants (Table 2), and these species can 
reduce soil moisture loss (soil evaporation and plant transpiration) due to their degree of coverage 
and drought-tolerance (Jiao et al., 2011; Hu et al., 2016). Finally, accumulation of plant biomass 
could offer food sources to soil animals (Table 1), and disturbance by animals could decrease BD 
and increase soil clay content during secondary succession (Hanél, 2003; Jing et al., 2014). Ren et 
al. (2017) also found that plant diversity (M, H, and E) synchronously increased with soil 
microbial diversity and biomass during restoration, and that microbial dynamics could alter soil 
texture (e.g., BD and clay). However, pH changed more slowly with successional time than SWC, 
BD, soil temperature, and clay (Table 3), which is consistent with the results of previous studies 
(An et al., 2008; Zhang et al., 2016). This suggests that secondary succession of abandoned 
farmlands has little effect on soil pH, since pH is the result of multiple synergies between biotic 
and abiotic factors (Criquet et al., 2000). 

In addition, our results showed that soil chemical properties (SOC, TN, TP, NH4-N, NO;-N, 
and AP) were coupled to changes in plant diversity, coverage, AB, and BB (Fig. 3a). Not 
surprisingly, the colonization of plants fixes CO» from atmosphere through photosynthesis, and 
returns C to soil in the form of litter and rhizodeposition starting from the initial period of 
succession (Ohtsuka et al., 2010; Xiao et al., 2021). This cause SOC to increase significantly with 
successional time (Fig. 2a). In contrast to SOC, measures of TN, NH4'-N, NOs-N, TP, and AP 
decreased significantly at 5 a plot compared with farmland (Fig. 2b—f), which is attributable to the 
cessation of input from fertilizers and absorption by plants (Zhang et al., 2016). As above- and 
below-ground biomass gradually increased, contents of TN, NH4'-N, and NO3-N increased 
significantly with successional time (Fig. 2b, e, and f) due to the increase in Leguminosae plants 
with N-fixing capacity causing an increase in N content of plant-soil system (O'dea et al., 2015). P 
is a more 'rock-derived' element and is fixed only slowly from atmosphere (Huang et al., 2013), 
and continuous increase in above-ground biomass of plants will absorb a large amount of P from 
soil (Chen et al., 2000), causing significant decreases in TP and AP contents with successional 
time (Fig. 2c and d). Although soil nutrients declined at the early stage of succession to some 
extent, they rose to farmland levels at 10 and 18 a plots and exceeded farmland levels at 30 a plot 
(Fig. 2), indicating that soil nutrient status can be improved with successional time (Wang et al., 
2009; Zhang et al., 2016). Net primary productivity of vegetation was enhanced during secondary 
succession (Table 1), causing increases in quantity of plant litter and dead roots, and ultimately 
improving soil nutrient levels (Peichl et al., 2012). Simultaneously, improving soil environment 
accelerated the decomposition of organic matter and then increased the accumulation of soil 
nutrients (Zhang et al., 2018b). Additionally, growth of plants controlled soil erosion, which also 
reduced losses of soil C, N, and P due to surface runoff (Zhu et al., 2021a; Luo et al., 2022). 

Development of plant coverage and biomass is the foundation of soil restoration (Ren et al., 
2017), and also strongly affects the dynamics of soil biological properties (Ren et al., 2016a). The 
results showed that soil enzymatic activities (SAC, URE, ALP, and CAT) and microbial 
biomasses (MBC and MBN) significantly increased with successional time (Figs. 4 and S3), and 
such alterations were also coupled to shifts in plant biomass and soil physical-chemical properties 
(Figs. 5 and 6). This indicated that secondary succession improves plant residues and soil 
microenvironment, and ultimately promotes soil biological activities. For example, increased 
plant residues during secondary succession contain a large amount of substrates that stimulate 
microbial growth and enzyme synthesis (Lucas-Borja et al., 2016). An improved soil environment 
(SWC, BD, and clay; Table 3) during secondary successional can promote the growth of 
microorganisms and contribute further to the development of microbial biomass and enzymatic 
activities (Yang et al., 2020; Guan et al., 2022). In addition, the development of below-ground 
biomass accelerates soil ventilation, promoting the release of exogenous enzymes by increasing 
rhizosphere exudates and ultimately altering microbial biomass and enzymatic activities (Zhang 
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et al., 2011; Zhang et al., 2018a). In contrast to microbial biomass and enzymatic activities, 
microbial respiration and qCO» increased significantly at 5 and 10 a plots compared with 
farmland but then decreased dramatically with successional time (Fig. 4g and h). As qCO» 
indicates the maintenance energy size of microbial communities and utilization efficiency of 
substrate, it is a sensitive index that reflects the influence of environmental factors on microbial 
activity (Liu et al., 2010). At early stage of succession, fast-growing vegetation requires a large 
quantity of nutrients, so microorganisms must accelerate fixation of nutrients to meet the 
vegetation's needs (Zhang et al., 2016). The results showed that qCO2 and microbial respiration 
increased when the rate of microbial turnover and energy consumption increased (Plaza et al., 
2004; Zhang et al., 2011). As succession progresses, the reduced environmental stress (lower BD 
and higher SWC) increases microbial efficiency and decreases qCO» because microorganisms 
need to invest little energy on maintenance (Ren et al., 2018; Zhao et al., 2018). This allows the 
soil to retain sufficient active organic matter to maintain good properties and sustainable 
utilization potential (Liu et al., 2010; Ren et al., 2018). 


4.3 Effects of plant community composition, species diversity, and dominant species on soil 
properties 


Plant community composition can reflect varying quality and quantity of plant residues (litter and 
rhizodeposition) provided to the decomposers (Yahdjian et al., 2017), and may also cause diverse 
micro-environmental conditions (Kardol and Wardle, 2010). Such resultant shifts are related to 
sources and turnover of soil nutrients and ultimately alter soil physical-chemical and biological 
properties (Xiao et al., 2017). The results showed that plant community composition was more 
closely related to soil physical-chemical and biological properties than species diversity and 
dominant species (Figs. 3 and 5), and indicated that plant functional group can better reflect soil 
properties. Li et al. (2018) found that plant functional groups representing a class of plant species 
with similar physiological and biochemical characteristics (N fixation and C fixation), which is 
more capable of dominating nutrient feedback and circulation of ecosystem than single species 
(Wu et al., 2017; Chou et al., 2018). Previous research has also shown that plant functional group 
composition can affect soil environmental conditions and nutrient supply to determine species 
diversity (Chen et al, 2021; Sidlauskaite et al., 2022), therefore, it can better reflect the 
interactive relationships between vegetation and soil than species diversity. In addition, the results 
showed that Leguminosae plants were the most important plant community factor associated with 
soil physical-chemical and biological properties (Fig. 6). Leguminosae plants play important role 
in promoting soil C and N accumulation and mineralization due to their capacity for N»-fixation 
by symbiotic root bacteria (Wu et al., 2017). Furthermore, Leguminosae plants contribute a net 
increase in litter N concentration due to their strong N fixation ability and low degree of N reuse 
from old leaves to new leaves (Li et al., 2012; Li et al., 2015). In addition, litter with higher N 
content contains more unstable complex, which is easily decomposed by microorganisms (Fanin 
and Bertrand, 2016; Zhang et al., 2018a; Zheng et al., 2021). Therefore, it can be concluded that 
Leguminosae plants promote the transfer of nutrients between plants and soil by accelerating the 
decomposition of litter, and increase the functional diversity and activity of soil microbial 
community by enhancing substrates inputs (Li et al., 2012; Hu et al., 2016). Increasing MBC and 
MBN and decreasing qCO» and microbial respiration over successional time confirmed the above 
deduction (Fig. 4). Fterich et al. (2014) also found that Leguminosae plants could increase soil 
organic matter and had positive effect on soil enzyme activity as well as C and N use efficiency of 
microorganisms. These variations in soil biological variables were significantly correlated with 
soil physical-chemical properties (Fig. 6), indicating that Leguminosae plants can change soil 
nutrients and environmental conditions by affecting soil microbial and enzyme activities. In 
contrast, co-improvement of soil physical-chemical and biological properties caused by legume 
species stimulates the development of mycorrhizal fungi, which strengthen plant nutrient 
absorption (Zhang et al., 2011), which lead to shifts in species diversity and plant community 
composition (Tables 1 and 2). Above all, plant functional composition can better reflect soil 
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physical-chemical and biological properties than species diversity and dominant species, with 
Leguminosae plant in particular being used to indicate the restoration status of vegetation and soil 
during secondary succession. 


5 Conclusions 


Plant communities transitioned from annual herbs to perennial species during 30 a of secondary 
succession. IV of Leguminosae and Gramineae significantly increased with successional time, 
whereas IV of Compositae significantly decreased. Moreover, plant community coverage, AB, 
and BB decreased at 10 a plot, and was then increased between 10 and 30 a plots, whereas species 
diversity significantly increased with successional time during first 18 a of secondary succession, 
then decreased at 30 a plot. Soil physical-chemical and biological properties (except TP and AP) 
significantly improved during secondary succession, even though some variables (e.g., TN, 
NH4-N, NO3-N, and SWC) degraded at the early stages compared with farmland. Changes of 
plant community composition had greater effects on soil physical-chemical properties and 
biological properties than species diversity and dominant species. Particularly, Leguminosae 
plants were the most important factor affecting soil physical-chemical and biological properties, 
and therefore can be used as indicators for observing the restoration status of vegetation and soil 
during secondary succession. This study provides suitable planting species for better vegetation 
restoration during secondary succession in the Chinese Loess Plateau. 
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Appendix 
Table S1 Geographical features and dominant species at different plots 
Plot Slope degree (°) Aspect (?) Elevation (m) Location Dominant species 
Farmland 26 North by east 45 1187.3 36?46'N, 109?15'E - 

5a 28 North by west 40 1152.4 36?45'N, 109?15'E A. capillaris, H. altaicus 

10a 30 North by east 55 1239.1 36°44'N, 109°16'E S. bungeana, A. capillaris 

18a 27 North by west 38 1195.6 36?45'N, 109?16'E L. dahurica, A. sacrorum 

30a 32 North by east 42 1251.2 36?44'N, 109°16'E A. giraldii, L. dahurica, B. ischaemum 


Note: A. capillaris, Artemisia capillaris Thunb.; H. altaicus, Heteropappus altaicus (Willd.) Novopokr; S. bungeana, Stipa bungeana 
Trin.; L. dahurica, Lespedeza dahurica Schindler; 4. sacrorum, Artemisia sacrorum Ledeb.; A. giraldii, Artemisia giraldii Pamp.; B. 
ischaemum, Bothriochloa ischaemum (Linnaeus) Keng. - indicates no dominant species. 


Type of enzyme 


Table S2 Methods for determination of soil enzymatic activities 


Detailed measurement method 


Soil catalase activity 


Soil catalase activity was determined by addition of 40 mL distilled water and 5 mL 0.3% H202 
to 2 g fresh soil. The mixture was shaken for 20 min (at 150 r/m) and filtered (Whatman 2V) 
immediately. Then the filtrate was titrated with 0.1 mol/L KMnO, under the conditions of 
sulfuric acid. Finally, the results were expressed as 0.1 mol KMnO,/(g-20 min). 


Soil saccharase activity 


Soil saccharase activity was determined using 8% glucose solution as substrates. About 5 g fresh 
soil was incubated with 15 mL substrates, 5 mL 0.2 M phosphate buffer (pH 5.5), and 5 drops of 
toluene for 24 h at 37.8?C. After incubation, the mixture was filtered (Whatman 2V) immediately 
and 1-mL aliquot was reacted with 3 mL 3, 5-dinitrylsalicylate in a volumetric flask, and then 
heated for 5 min. Soil solution in the flask was quantified in an ultraviolet spectrometer 
subsystem (UVS) at 508 nm when it reached room temperature. Finally, results were also 
expressed as mg glucose/(g-24 h). 


Soil urease activity 


Soil urease activities was routinely determined using 10% urea solution as substrates. About 5 g 
fresh soil was incubated for 24 h at 37.8°C with 5 mL citrate solution at pH 6.7 and 5 mL 
substrates. The reaction mixture was then diluted to 50 mL with distilled water. After incubation, 
the mixture was immediately filtered and 1 mL supernatant was treated with 4 mL sodium phenol 
solution and 3 mL 0.9% sodium hypochlorite solution. The released ammonium released from 
urea hydrolysis was quantified in an ultraviolet spectrometer subsystem (UVS) at 578 nm. 
Results were expressed as mg NH4'-N/(g:24 h). 


Soil alkaline phosphatase 


activity 


Fig. S1 


successional time 


Soil alkaline phosphatase activity was determined by addition of 10 g fresh soil, 2 mL toluene, 
10 mL disodium phenyl phosphate solution, and 10 mL 0.05 M borate buffer. The reaction 
mixture was incubated for 2 h at 37.8?C. After incubation, the mixture was immediately filtered, 
then the filtrate was treated with 0.5 mL of 2% 4-aminoantipyrine and 8% potassium 
ferrocyanide; the phenol released was determined in an ultraviolet spectrometer subsystem 
(UVS) at 510 nm. Results were expressed as mg phenol/(g-24 h). 
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Fig. S2  Nonmetric multidimensional scaling (NMDS) analysis of soil physical-chemical properties over 
successional time 
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Fig. S3 Nonmetric multidimensional scaling (NMDS) analysis of soil biological properties over successional 
time 


